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We study the rate of quantum phase slips in an ultranarrow superconducting nanowire ex¬ 
posed to weak electromagnetic radiations. The superconductor is in the dirty limit close to the 
superconducting-insulating transition, where fluxoids move in strong dissipation. We use a semi- 
classical approach and show that external radiation stimulates a significant enhancement in the 
probability of quantum phase slips. This can help to outline a new type of detector for microwave 
to submillimetre radiations based on stimulated quantum phase slip phenomenon. 


I. INTRODUCTION 

Quantum phase slip junctions are exact dual counter¬ 
part of the Josephson junctions. Recently these junctions 
have been successfully realized in ultranarrow supercon¬ 
ducting nanowires, where quantum phase slip replaces 
tunneling Cooper pairs mm- These nanowires are non¬ 
linear elements performing similar physics as Josephson 
junctions with the roles of superconducting phase ip and 
charge q being interchanged mm- This duality has been 
the motivation behind many of the recent applications 
of quantum phase slip (QPS) elements [5]-[8]. These ele¬ 
ments have found interesting implications for fundamen¬ 
tal metrology and information technology, for instance as 
photon pulse detectors, quantum current standard, and 
quantum bits [IHIQ]. 

In a superconducting nanowire with small cross sec¬ 
tion, the supercurrent is determined by the phase differ¬ 
ence (p between two ends of the nanowire from the saw¬ 
tooth relation Ig = ^opf27TL with L being nanowire ki¬ 
netic inductance and ^0 = h/2e. In temperatures much 
lower than the superconducting critical temperature (i.e. 
T <^Tc) quantum fluctuations may suppress the modu¬ 
lus of the order parameter in a region and turn it from 
superconductor to normal metal. This enables the super¬ 
conducting phase to slip by 2n7r, with integer n, without 
any energy compensation. An individual phase slip takes 
place in a normal core, similar to the normal core of a 
magnetic flux vortex, therefore we can assume the core 
size is roughly the coherence length ^ mm- QPS event 
takes place for a short period of time that is maximally of 
the order of inverse of superconducting gap h/2A. Simi¬ 
lar effect happens close to Tc due to thermal fluctuations 
of the order parameter [13] . 

In superconducting nanowire made of clean materials 
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with low normal resistance R, quantum phase slips rarely 
take place. To enhance the slip rate a nanowire should 
be made of highly disordered amorphous superconduc¬ 
tor, which is in the dirty limit, with large R [5j[7]. There 
is not a well-understood theory to describe the super¬ 
conductivity in near superconductor-insulator transition 
(SIT). A candidate theory [T4l|T5] proposes superconduc¬ 
tivity at high disorder is maintained by a fragile coher¬ 
ence between electron pairs, which is characterized by an 
anomalous binding energy. If pairs are localized, they en¬ 
ter an insulating state, and if condense, a coherent zero- 
resistance state emerges. Based on this theory supercon¬ 
ductor in SIT have regions of localized BCS-condensates 
nearly separated in different lakes [16] . The cores of QPS 
can coherently tunnel across superconducting regions and 
avoid dissipation. This is similar to the Cooper pairs that 
tunnel across a Josephson junction without much dissi¬ 
pation [miiH]. 

The voltage across the nanowire is known to be pe¬ 
riodic in charge of the crossing fluxoid; i.e. V = 
Vb sin(27rg/2e). Individual phase slips in nanowires can 
be observed when a large bias voltage is applied on the 
wire. Under such bias voltage, effective potential be¬ 
comes a tilted washboard with more slanted slop in larger 
bias. Depending on temperature, there are two general 
scenarios for the dynamics of a fluxoid. Close to the crit¬ 
ical temperature Tc, fluxoid particle gains energy from 
thermal activation and overcomes potential barrier to slip 
across the wire m- Quite differently, in low temperature 
T Tc fluxoid particle becomes frozen in a minima of 
the washboard potential. The minima are called ‘zero- 
current states’ where Coulomb blockade occurs [ME]. 
Vacuum fluctuations help the particle to tunnel into the 
barrier and slips away. Coherent tunneling is possible 
between two zero-current states where a quantum vari¬ 
able (phase or charge) has minimum fluctuations. Such 
coherent tunnelings have been previously observed in the 
superconducting-insulating transition limit [T9]-[2T]. 

Exposing nanowire to strong electromagnetic radia- 
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tion produces Shapiro steps [22l [23] in the current- 
voltage dependence, which has been observed in reference 
[24|. However, in many applications of superconducting 
nanowires, such as in qubits and photon detectors, weak 
radiation is applied where phase locking cannot occur. 

In this paper, we qualitatively study the effect of a 
weak alternating electromagnetic field on the quantum 
phase slip rate in ultranarrow superconducting nanowire, 
where the width of the nanowire is smaller than the su¬ 
perconducting coherence length, i.e. r < ^. We consider 
the nanowire is in the insulating phase. Our method 
is to map this problem into its well-studied analogue 
in Josephson junction in proper regime. We use the 
semiclassical quantum mechanical approach developed 
by Ivlev and Mel’nikov [25II27] in studying quantum tun¬ 
neling in a high-frequency field to our problem. Similar 
to a Josephson junction under weak time-harmonic ra¬ 
diation EZlllH], we expect a significant enhancement in 
the stimulated phase slips at zero temperature. We show 
that a fluxoid gains energy from radiation and tunnel 
into the barrier more often than usual and slips away. 
This leads to the super-exponential enhancement in the 
rate of such ‘stimulated quantum phase slips’ (SQPs). In 
certain nanowires, this can result in larger DC resistivity 
with minimal fluctuations in a dynamical variable. 

The enhanced escape from the zero-current state stim¬ 
ulated by weak irradiation has significant practical im¬ 
portance. Since the observation of zero-current state re¬ 
quires a high impedance environment, our study in here 
will be confined to the highly dissipative cases. Our 
model suggests that the QPS in certain nanowires at low 
temperature can significantly be amplified in microwave 
to THz radiation. The feasibility of this detector at the 
typical frequency of 0.3 THz using conventional materials 
will be discussed. 


II. STIMULATED QUANTUM PHASE SLIPS 

A superconducting nanowire with QPS is the dual to 
a Josephson junction with charge and phase (as well as 
current and voltage) interchanged. In Josephson junc¬ 
tion, a Cooper pair tunneling across the junction picks 
up a phase exp(±i(^(t)) corresponding to the supercon¬ 
ducting phase This induces a coupling energy E = 

Ej{l — cos ip). The current is defined I = (2elK)dE jdip. 

Analogously for a nanowire similar relations can be 
derived. A QPS fluxoid picks up a charge phase when 
tunneling exp(±i(5) with Q = 2'Kql2e being a dimen¬ 
sionless charge parameter. Therefore the QPS energy 
becomes ^ — cos Q{t)). The voltage is defined 

V = dE/dq. Phase slips may take place everywhere in 
the wire whose induced current depends on the wire in¬ 
ductance. Therefore a narrow superconducting nanowire 
can be modelled as a voltage in series with an inductance 
as shown in the left part of Fig. In the figure, the dis¬ 
sipation is modelled by a resistor and the AC and DC 
bias voltages are sources in series with the wire. This 



FIG. 1. The schematic circuit of a QPS junction includ¬ 
ing an ideal QPS element (a superconducting nanowire), the 
dissipative element R, the bias voltage Vdc and the driving 
source Vac- L denotes the dominant kinetic inductance. 

circuit is built based on laia. The inductor L is the 
total of the kinetic inductance {Lk) and the geometric 
inductance (Eg) of the circuit. Since, in superconducting 
nanowire, the kinetic inductance is much larger than the 
geometric inductance, we have L ^ Lj^. In the circuit of 
Fig. [^voltage is V = Vo sin {27Tqf2e) + Lq J- Rq with the 
QPS and inductance energies 

Es = 2ero/27r, El = ^1/21^. (1) 

where Vo is the voltage scale of QPS energy Es- 

In the nanowire, a crossover from insulator to a su¬ 
perconducting inductor takes place when the inductance 
energy El is increased beyond QPS energy Es- In 
the superconducting phase El ^ Es^ the fluxoid en¬ 
ergy E = Es{l — cosQ) + EL{(t)f) is dominated by the 
parabolic inductance energy EL{(t)f) associated with in¬ 
duced phase (pf. The parabola associated with differ¬ 
ent winding integer n cross at certain energies where the 
small energy of QPS provide an avoided crossing gap. 
This makes nanowire energy to be multivalued in sep¬ 
arated energy bands, similar to a capacitive Josephson 
junction. 


In the opposite regime where Es ^ El the wire en¬ 
ergy is dominated by QPS energy — cosQ) which 
oscillates in charge. Consider that the charge under¬ 
goes a fluctuation around its macroscopic value q so that 
the stochastic charge is q = q ^ 5q. Following the ana¬ 
logue discussion for a Josephson junction phase, (see [29F 
[3T] and references therein) the charge fluctuation cor¬ 
responds to the effective current fluctuations across the 
wire 

= [ SI{t')dt'. (2) 

Jo 

A difference between QPS and thermally activated 
phase slips (TAPS) is that the dissipation in TAPS is 
due to stochastic energy activation in high tempera¬ 
ture while QPS allows tunneling between distinct zero- 
current states. The latter is similar to zero-voltage states 
in Josephson junction [T9]-[2T]. For a nanowire in the 
insulating phase the quantum tunneling between zero- 
current states takes place without current fluctuations 
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Therefore Eq. <§ shows that the charge in fluxoid 
behaves semiclassically, whereas superconducting phase 
can be subject to large fluctuations [4]. 

The semiclassical charge associated with QPS fluxoid 
in the nanowire depicted in Fig. [Revolves in the following 
way: 

(cos Q -ko -ki cos fit) = 0, (3) 

with O is the frequency of the driving voltage and 

7] = R/L, Up = y27rVb/2eL, ki = V^/Vo- (4) 


with the index i = 0 (1) corresponds to DC (AC) volt¬ 
age. The definition of the plasma frequency ujp is com¬ 
patible with the definition based on the duality: hujp = 
^/2 EsEl^ which is similar to the equation of RCSJ model 
of Josephson junction [32l |33] with high-frequency driv¬ 
ing field. For simplicity in writing Eq. (§ we had shifted 
Q Q IT 12 in order to have applied and QPS voltages 
in phase. 

Our aim is to study the possibility of utilizing nanowire 
as a detector for time-harmonic radiations, therefore we 
restrict ourselves to weak alternating fields, i.e. ki 1. 
This makes our problem to be different from the physics 
of the Shapiro steps [22j where the phase of nanowires (or 
its dual Josephson junction) is locked to the frequency of 
the driving field frequency and constant voltage steps are 
observed. For weak time-harmonic fields the driving force 
is very small and the wire is in zero-current state with 
ko < 1. The most important result is that the smallness 
of ki does not necessary mean that the its effect on the 
charge dynamics is small. In fact as we will show a weak 
time-harmonic field can significantly affect the wire by 
increasing the rate of QPS (see also Appendix [A|). 

In the limit of weak dissipation, the tunneling rate of 
fluxoid can be studied using ID quantum mechanics of 
the Lagrangian associated to Eq. (§ subject to R = 0. 
However we are interested to study the decay of the zero- 
current state in the limit of strong dissipation because 
resistance is large in wires with QPS effects. A dual ef¬ 
fective theory has been developed for dissipative coherent 
tunneling in Josephson junction by Caldeira and Leggett 
pn [TS] , and Larkin and Ovchinnikov [34] . In those the¬ 
ories, dissipation is modelled as the coupling to bosonic 
degrees of freedom. The low-energy effective action of 
the system is derived to properly take account of the dis¬ 
sipation. 

The semiclassical theory of Josephson junction exposed 
to weak alternating current in EH El guides us to study 
the charge dynamics of radiation-stimulated quantum 
phase slips (SQPS) in nanowires with strong dissipation. 
Effectively the evolution of semi-classical charge that tun¬ 


nel across the wire is: 


cPQ ^ 2, n 

-^+Up{cosQ 


ko — ki cos Vtt) 


—2iiiri 


¥)( 


sin[(g(^) - Q{ti))/2] ^ ^ 
c ^ {TTkBT{ti — t)/h) 


(5) 


where the contour C is shown in Fig. |A.8| and the princi¬ 
ple value of the integral is implied (for general discussion 
of the method, see Appendix 0- 

In the limit of our interest, the nanowire is strongly 
dissipative rj ^ ujp. For semi-classical description to be 
valid, it is required that Es ^ hCl El El. Also the 
applied DC voltage is close to Vb, be. Vo — Vdc ^ Vb, 
therefore the term with second derivative in Eq. (§ can 
be omitted. Also, we can assume that the exchange of 
energy between the wire and its environment takes place 
in the shortest time, thus the argument of sinh in the de¬ 
nominator of Eq. <§ can be replaced by its lowest order 
sinhx ~ X. Technical analysis of this integration over the 
contour C shows that the integral tends to zero except 
that at the singularity t = ti where its proper residue 
must be counted (see Eq. (18) in Ref. [27]). Therefore, 
in the lack of alternating radiation Eq.® in the regime of 
interest effectively reads: —r]dQ/dt + (cosQ — 1) = 0, 
which has the following solution 


Q(t) = iln 


t - ITs 
tElTs' 


Ts 


T] 




(6) 


with Ts being the time of under-barrier motion. 

In a system described by the classical action S = 
—i Jfj Cdt with C being Lagrangian, the probability of 
quasiclassical tunneling is T = exp (—S'). Regarding the 
alternating voltage being a small perturbation ^ 1, 
we can rewrite the action in the form of S = So + Si with 
So being the action in the lack of alternating field and Si 
is linear in Vac, [35] . 

Let us assume the QPS probability in a nanowire with 
strong dissipation and DC voltage Vdc about Vo is de¬ 
noted as To- Above results easily show that in the 
presence of a weak high frequency radiation hitting the 
nanowire, the probability of QPS in the wire will change 
from To to F in the following form: 


r(F4C'A) = To exp 


4eVAC 

hfl 


sinh(firs) , 


(7) 


Eq.Q is the main result in this paper. It indicates 
that an alternating driving field with certain frequency 
and voltage can trigger occurance of a large number of 
QPS’s in a proper nanowire at low temperature. For 
instance if a wire with dissipation factor rj/ujp = 10 is 
driven by a weak time-harmonic radiation of the real- 
tive amplitude AeVAC/^p = 5 x 10“^ and frequency 
D = cjp, the QPS rate increases by a factor of about 250 
times. A more careful analysis shows that this result is 
valid for nanowire temperature T < To with To being 
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the crossover temperature bet ween quant um and ther¬ 
mal activation regimes Tq = y^( 1 — /co)/2 jirkBr], 

m- 

From Eq.Q one can see that the bigger the normal 
resistance R is, the larger the rate of SQPS becomes. In¬ 
tuitively this can be understood from the definition of 
underbarrier time in Eq. (H- Upon increasing dis¬ 
sipation the under-barrier time grows larger. According 
to the semiclassical description of quantum tunneling, 
see Appendix during quantum tunneling time pa¬ 
rameter becomes imaginary. This changes the bounded 
function of alternating potential in eq. (§ into the un¬ 
bounded function cosh Ur. The longer a fluxoid stays 
under the barrier, the more energy it absorbs from the 
alternating potential and this causes stimulation of quan¬ 
tum phase slips. 

The QPS rate in the lack of an time-harmonic drive 
has been estimated by Mooij and Harmans in [6] to be 
nearly To ~ Es/^- This in addition to substituting Eqs. 
0 and 0 simplifies Eq. 0 into: 


r(14c,^^) = Vo/^o)exp 


4:eVAc 

hfi 


sinh ( 



(8) 


Given that the QPS rate in the absence of time-harmonic 
radiation is proportional to Vq, one of the features of 
Eq. 0 is that the super-exponential enhancement of 
QPS rate is inversely proportional to the rate Tq, thus for 
small rate Tq the exponential enhancement of QPS in the 
presence of high-frequency field is more significant. This 
enhancement is only due to the stimulated excitation in 
the zero-current states in the presence of external drive. 

In the absence of the Udc, there is no tilt in potential 
and the rate of fluxoid crossing to right or left are equal. 
As a result the average current becomes zero (7 = 0). 
However, in the presence of a positive value for DC volt¬ 
age bias, the average current is given by: 


/ = 2e(r^-r^) 


(9) 


with F^ (r^) the rate of crossing to the right (left) 
where the potential barrier decreases (increases). In the 
case the bias voltage is close to the critical voltage (Vb), 
F^ the crossing F^ dominantly exceeds that of the op¬ 
posite direction. Hence F = F^ and 


/ = 2er(U4c,^), 


( 10 ) 


where T{Vac,^) is given by Eq. 0 According to Eq. 
(10), the influence of high-frequency weak irradiation 
on superconducting nanowire biasd by the DC voltage 
Vq — Vdc ^ Vb is observable by measuring the crossing 
current. 

The quality factor in nanowire Qs- is defined as: 


Qs = ‘^. ( 11 ) 

v 

In low quality factor QPS nanowire the dissipation is 
strong. The larger r] leads to longer under-barrier time 


and consequently the enhancement of SQPS rate expo¬ 
nentially increases. The increase in the under-barrier mo¬ 
tion due to higher dissipation can be seen in Eig. 1^. 
Therefore, we expect that a low-Q^ nanowire to be a 
better candidate for observing tunneling enhancement. 

A comment on the range of validity of the method 
we used in this section is in order. As it is seen from 
Eq. 0, the enhancement in the tunneling probabil¬ 
ity for fir] ^ ujp is itself an exponentially large factor 
(^ eU 4 C'(^^)~^ exp(U77/cc;^)). This indicates that the 
range of the validity of the semi-classical approach in 
this case is limited to cVac r\j (HQ) exj[){—Qr]/ujp). Be¬ 
yond this, higher order correction in terms of Vac to the 
calculations is required [36] . 

An alternative approach in studying the QPS rate in 
superconducting nanowires under high-frequency radia¬ 
tion would be to use the the effective action method 
developed by Golubev and Zaikin [TJ [37H39] in a non¬ 
equilibrium setting. There are some challenges associ¬ 
ated with this approach that are studied in [40] . 


III. PROPOSAL FOR QPS-BASED 
ENERGY-RESOLVING HIGH-FREQUENCY 
RADIATION DETECTOR 

The exponential enhancement of the probability of the 
quantum tunneling observed in Eq. 0 can be exploited 
in designing detectors of microwave to THz radiation that 
are capable of determining the frequency of the incom¬ 
ing energy. In this section we introduce a new type of 
high-frequency detectors based on enhancement of QPS 
phenomenon in superconducting nanowires. In this pa¬ 
per, the working principle of this type of detector will be 
discussed and much important engineering details like 
impedance matching will not be addressed. 

Our proposed detector is made of a low quality fac¬ 
tor, i.e. low Qs, QPS junction that is voltage biased 
close to the critical voltage. An antenna is the source 
of the high-frequency voltage and is placed right across 
the superconducting nanowire. The current in the loop 
is measured constantly, the change in the current and 
the amplitude determines the presence of the detected 
radiation. A schematic of this system is shown in Eig. 

0 . 

The superconducting nanowire is placed in the gap be¬ 
tween two parts of the antenna. This will guarantee that 
the maximum Vac is induced along the nanowire. Other 
elements like the resistance and the voltage-bias source 
are placed outside of the antenna in the loop. 

Presence of the radiation results in the decay of the 
zero-current state of the QPS junction which causes a 
change in the current of the circuit. Depending on the 
design parameters, the detection of the change in the cur¬ 
rent might be hard to achieve. A lock-in amplifier or a 
SQUID can be used for current monitoring in case the 
current change is difficult to be monitored with conven¬ 
tional methods. 
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FIG. 2. The schematic of a QPS high-frequency detector is 
shown. The red segment in the middle, is the superconducting 
nanowire. A broadband bow-tie antenna collects the high- 
frequency field. The resistance R adds dissipation to lower 
the quality factor of the QPS junction. 


The resistance R plays the important role of reduc¬ 
ing quality factor {Q) of the QPS junction. Its value is 
chosen such that the required enhancement in Eq. 0 
is achieved, which depends on other parameters of the 
system. 


A. Design Parameters 


As an example, in this section, we investigate design 
parameters for a Q/27r = 0.3 THz detector. We will study 
different superconducting materials to evaluate the appli¬ 
cability of them in our design. Based on these properties, 
parameters of the detector can be estimated. 

We assume that the QPS energy is related to the QPS 
rate according to Es = hTqps [6]. The QPS rate from 
the Golubev-Zaikin theory misMi is given by 

^ AR„X^ ( RnX\ 

= ( 12 ) 


where R^ is the normal resistance per unit length of the 
superconducting nanowires. The two constants Ci and C 2 
account for uncertainties in derivation of Eq. (12) which 
are of order one. We set ci = C 2 = 1. Although Eq. (12) 
is given by Golubev-Zaikin theory, the factor 0.3 in the 
exponent is adopted from the fit of experimental data to 
the Giordano model in the work of mm- 

In order to choose the appropriate material and pa¬ 
rameters for the detector, we study properties of four 
different materials NbSi, InO^^, NbN and Ti. Properties 
of these materials are listed in Table HI The coherence 
length ^ for superconducting nanowire is related to the 
bulk parameter through 


^ ~ 0.85^Cbulk^o, (13) 

where, Iq is the mean free path of the electrons, narrow 
superconducting nanowires are always in the dirty limit. 


Eig|^ shows Es for four different materials as a func¬ 
tion of normal resistance per length. The resistance per 


TABLE I. Material properties of the superconducting 
nanowires used for simulations. 


Material A 

[meV] ^ 

[nm] 

NbSi 

0.18 

15 

InOx 

0.41 

20 

NbN 

1.6 

4 

Ti 

0.06 

80 


The data for NbSi, InOx, NbN and Ti are adopted from m, m, 
m and [2S] respectively. 



FIG. 3. The QPS energy as a function of the normal resis¬ 
tance per length for four different materials NbSi, InOx, NbN 
and Ti is shown. The length of the nanowire is A = 2/im. 
Parameters are listed in Table [U 


length Rn determines the cross-section area of the super¬ 
conducting nanowire. Since Rn is inversely proportional 
to the cross-section area of the nanowire, the higher Rn 
indicates thinner nanowires that support smaller current 
which makes the operation more difficult. 

Another important energy scale in QPS junctions is the 
kinetic inductive energy Ep which plays an important 
role in the dynamics. The kinetic inductive energy is 
given by 


F 

^ 2Lfe’ 

where the kinetic inductance is found from 

hRjsf 


Lk = 


ttA 


(14) 


(15) 


In Eq. (15), Rn is the total normal state resistance of the 


superconducting nanowire which is given hy Rn = XRn, 
where X is the length of the superconducting nanowire. 
In Eq. (14), the geometric inductance and external in¬ 


ductance are assumed to be much smaller than the kine¬ 
matic inductance Lk of the superconducting nanowire. 
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FIG. 4. The inductive kinetic energy as a function of the 
normal resistance per length for four different materials NbSi, 
InOcc, NbN and Ti is shown. The length of the nanowire is 
X = 2/im. Parameters are listed in Table [I] 


FIG. 6. The ratio of EsIEl as a function of the normal 
resistance per length for four different materials is shown. 
For the charge number to be the good quantum number, it is 
required that Es > 4 :El- 



FIG. 5. The plasma frequency as a function of the 

normal resistance per length for four different materials NbSi, 
InOcc, NbN and Ti is shown. The Rn determines the dimen¬ 
sions of the superconducting nano wire. The length of the 
nanowire is X = 2/im. Parameters are listed in Table 


In Fig. [^and Fig. [^the inductive kinetic El energy 
and the plasma frequency uOp of the four nanowire as a 
function of normal resistance per length are shown. 

Since the superconducting nanowire is intended to be 
working in the regime where charge is a good quantum 
number; this requires that at least Es > ^El- The ratio 
of EsIEl is shown in Fig. The acceptable region of 
parameters is anywhere above 4. 

Arbitrarily, we choose the length of the superconduct¬ 
ing nanowire to be X = 2/im. To satisfy the conditions 


of the semiclassical approach of the previous section, i.e., 
Es '> hft, we choose Es = 3 THz for the operating fre¬ 
quency of ft/27r = 0.3 THz. This leads to the critical 
voltage of Vq = ^Es = 39mV. 

Assuming the induced alternating voltage collected by 
the antenna has the amplitude of Vac = 100 nV, then 
AcVac/^^ ^ 10“^. Therefore, from Eq. (|^, to have a 
significant enhancement it is necessary to have 

^ > sinh“^( 10 ^) ^ 8 , (16) 

ujI 


where the dissipation 77 with dimension radian per second 
is defined in Eq. Q. Fig. I shows that InO^^ could be a 
possible candidate. In order to obtain Es = 3 THz, from 
Eigs. iilH and we obtain the following parameters 
for InOx nanowire: Rn ~ 0.025 K Ohm/nm, El ~ 120 
GHz, ujpl27T ^ 0.84 THz, and EsIEl ~ 24. Using these 
parameters in Eq. (16), in order to have significant en¬ 
hancement the total resistance of the circuit needs to be 
much larger than 490 K Ohm. Eor an on chip resistance, 
NiCr thin-film resistors can be used. 

One of the significant advantages of the proposed de¬ 
tector would be the simplicity of the fabrication using 
one-layer lithography. The superconducting nanowire 
and the antenna can be fabricated lithographically on 
the same substrate. Avoiding shunting parasitic capaci¬ 
tances might be challenging that requires extra attention. 

The width of the QPS element is in the order of 10 
nm to 20 nm, therefore, a large number of them can be 
fabricated in parallel which makes them good candidates 
for applications that require many elements like imaging 
or for higher detection and coupling efficiency. 

In the proposed method the detector is made of narrow 
nanowires with width smaller than the coherence length, 
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therefore the presence of vortices can be ignored. This 
is because structures smaller than 4.4^ can not support 
vortices. The absence of vortices might improve the noise 
performance of the proposed detector. 


IV. CONCLUDING REMARKS 

We studied the stimulation effect of a weak high- 
frequency field on the zero-current state tunneling of flux- 
oid particle. The approach chosen was to use the duality 
transformation between Josephson junction and a QPS 
junction to map the dynamics of QPS charge in a circuit 
model. Then we studied the effect of high-frequency al¬ 
ternating field on the coherent tunneling rate. The sim¬ 
ilar problem has been studied for the case of Joseph¬ 
son junction using semiclassical physics [25J[27] which 
we adopted for the case of QPS junction. We observed 
that in a strongly dissipative superconducting ultranar¬ 
row nanowire, a high frequency field can enhance the 
probability of quantum tunneling super-exponentially. 
Interestingly we find that the enhancement of SQPS rate 
is more pronounced in wires with small non-stimulated 
QPS. This result will help to predict that quantum phase 
slip qubits should be better-working in the presence of 
weak driving field. 

The rate enhancement and its driving frequency de¬ 
pendence can be exploited in designing energy-resolving 
high-frequency detectors. We outlined a new type of 
high-frequency detectors based on the QPS phenomenon 
in superconducting nanowires. The basic physics and 
design for such a detector was introduced. We then in¬ 
vestigated the possibility of such realization using the 
materials used in studying of QPS in superconducting 
nanowires. It was shown that the theoretical restrictions 
can be met by choosing correct design parameters. 
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Appendix A: Brief Review of the Quantum 
Tunneling in Time-dependent Potentials 

In this appendix, we review the effect of a high- 
frequency field on quantum tunneling in the semi¬ 
classical descript ion [36]. The approach will be based on 
the method developed in [25ll27] . Eirst, a brief introduc¬ 
tion is given to the semi-classical approach to quantum 
mechanics and quantum tunneling. 

The semi-classical description is obtained from the sta¬ 
tionary path approximation of the Eeynman path integral 


approach to quantum mechanics. The stationary path 
of a Eeynman path integral which is obtained from the 
variation of the action yields the Newtonian equation of 
motion (EOM): 


SS = 0 


EOM 


(Al) 


This relation is familiar in classical mechanics for ener¬ 
getically allowed region; however, the natural question 
that arises is that: is this method applicable to energet¬ 
ically forbidden regions like in quantum tunneling? The 
answer is “yes”; however, it requires allowing the time to 
acquire an imaginary part 02]. To see this, let’s consider 
the action of a point particle with mass m in a potential 
V{x). The action can be written as: 


S = 




^ hr -V{x) + E^ 


2 \dt 

The equation of motion is found to be 
d^x dV(x) ^ 


dx 


and the total energy is given by 


(A2) 


(A3) 


E=^ + Vix) 


(A4) 


where the mom entu m is defined as p{t) = mdx/dt. By 
integrating Eq. (A4), the required time t for the particle 
to reach infinity from point x is given by 


t{x) 


f 


dx'\/m 


^2{E-V{x')) 


(A5) 


Erom Eq. (A5), it is seen that as long as E > V{x) the 
time remains real, but ioi E <V (x) it acquires an imag¬ 
inary part and becomes complex. Therefore, by allowing 
complex time, classically forbidden regions can be stud¬ 
ied in the semi-classical approach. The tunneling of a 
point particle with energy E coming from left from the 
potential V{x) is shown in Eig. A.7 Eor x > X 2 time is 
real, because E > V{x); however, for xi < x < X 2 , the 
time goes in the imaginary direction. Eor x < xi time 
becomes complex, t + iro, where the constant imaginary 
part is 


ro = 


f 

J X\ 


dx^/m 


^2{V{x)-E) 


(A6) 


Therefore, for a tunneling path of a particle moving 
from left to right, the time evolution is depicted in Eig. 
|A.8| by contour 

According to the semi-classical description, the tun- 
nelin g amplitude is found by calculating the action Eq. 
(A2) along contour (7+. In order to find the tunneling 


probability amplitude the contour C- needs to be added. 













(A9) 



FIG. A.7. Potential barrier for a particle moving from left 
to right with energy E. Classical t urni ng points are indicated 
by xi and X 2 - According to Eq. (A5), for xi < x < X 2 , the 
time becomes complex. 




to + iTo 
- >1 


c_ 


-> t 


^0 “ 

FIG. A.8. The integration contour for the quantum tun¬ 
neling probability. The vertical sections correspond to the 
underbarrier motion. 


where the property x{E) = x*{t) has been used. There¬ 
fore, the tunneling probability with exponential accuracy 
is given by 


r « exp {-So), 

So = —i / dt 

Jc-+c+ 


m f dx ^ -r./ N ^ 

2 U' 


(A7) 


where x(t) is th e so lution to the classical equation of 
motion, i.e. Eq. (A3), along the contour. 


The horizontal segments of and C- cancel each 
other and only vertical segments corresponding to the 
under barrier m otion survive. Using Eq. (A3) in the 
exponent of Eq. (A7) we get 


: exp 


= exp I —2 


/ r-iTQ \ / 

li dt mx^ j = exp I 2im 

J dx\/2m{V{x) — E)^ , 


dxx 


(AS) 


which is the well-known WKB result in quantum me¬ 
chanics m- 

If the system is in thermodynamic equilibrium before 
the tunneling starts, then the tunneling probability needs 


to be statistically averaged over E 

E 


(r) = / 


dEexp 


knT 


-S{E) 


where E is given by Eq. (A7). The largest probabil¬ 


ity of tunneling occurs for energies that minimizes the 
exponent in Eq. (|A9[) and is given by 


dS{E) 


1 


(AlO) 


dE ksT 
The action for the underbarrier motion is given by 

rX2 

S{E) = 2 dx^/2m{V{x) - E), (All) 

J Xi 

and the energy derivative of the action yields 
dS{E) 


dE 


= -2 


/ 

xi 


dx 


^2{V{x) - E) 


= -2to (A12) 


where tq is the time of the under barrier motion given 
by Eq. (A6). Comparing Eqs. (AlO) and (A12) reveals 
that 


1 


^0 = 


2kBT 


(A13) 


Eq. (A12) determines the energy of the tunneling par¬ 


ticle. Therefore, in equilibrium, the proba bility of tun¬ 
neling through the barrier is given by Eq. (A7) for real 
trajectories that satisfy Eq. (A13). The real trajecto¬ 


ries condition comes from the analysis that shows that 
the time-averaged probability for semiclassical processes 
is entirely determined by real trajectories [25] . 

The semi-classical method in which time can take 
on complex values is suitable for generalization to in¬ 
clude tunneling from time dependent potentials, tun¬ 
neling from periodically modulated potential barriers is 
the most common application of this method and since 
in this paper we are interested in sinusoidal alternating 
field we restrict this section of this review to potentials 
of the form [T8l[25H27] 


U{x, t) = V{x) + £x cos Q,t. 


(A14) 


According to the semiclassical description, the linear in 
the field-strength correction to the tunneling probability 
is then given by 


Si = —i£ / dt x{t) cos 

Jc-+c^ 


(A15) 


where x{t) is the solution of the unperturbed equation 
of motion Eq.(A3) along the contour. During the under¬ 
barrier motion in the vertical segment of Fig. ITU time 
is imaginary t = ir and therefore the equation of motion 
becomes 


d‘^x dV(x) 
m—^ -^ = 0. 


dr^ 


dx 


(A16) 
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which in comparison to Eq. (A3) can be interpreted as 


the classical equation of motio n in the inverted potential. 
The contour in Eq. (A15) may be shifted to entails 


the singularities of the integrand. This enables calculat¬ 
ing the integral based on the singularities of the x {r). 
Therefore, the general trend of the Ai in Eq. (A15) de¬ 


pends on the specific form of the potential. In some cases 
Si can be exponentially large, which is the case we en¬ 
countered in this study. In order to have exponential 
enhancement, it is necessary for the function defined as 


h{x) = y/E - v{x), 


(A17) 


to have singularities off the real axis in the x plane 
Assuming V(x) has singularities of the form 


V(x) 


k(x — Xs)^, a < 0, X 
Kx^, a > 0, X ^ oc. 


(A18) 


then, the solution to Eq. (A3) near x^ is of the form 

x(t) =x,+ [- A(2 - - tsf] ^ > (A19) 

L 2m J 

where ts is the complex time that takes going from X 2 in 
Fig-iyg to Xg given by 


= 


f 


dx* 


^2{E-V{x'))' 


(A20) 


By comparing to Eq. (A6), the time of the underbarrier 


motion tq has the same order of the magnitude as 

rs=lmts. (A21) 

but for analytical potentials always rs < tq. Eor Qxg ^ 
1, the main contribution to the integral in Eq. (A15) 


comes from branch-cut section in the vicinity of the sin¬ 
gular points Ts and r*. Therefore, the transition proba¬ 
bility is given by 


r(f, t) = To exp (ai cos{Qt )), 


where 


ai = 


2tt£ 

rJ ^ ] 

-1 

do 

1 

to 

IT 

U-2; 


2mQ? 


(A22) 


exp(f]rs). 


(A23) 

where T^; is the Euler Gamma function and Tq is the 
tunneling rate in the absence of the alternating field. The 
time averaging of Eq. (A22) gives 


r(f) = n 


1 


y27rai 


^p(ai). 


(A24) 


Eq. (A24) shows that the semiclassical description is 
relevant only when 


5'o > ai > 1. (A25) 



FIG. A.9. Dissipative quantum tunneling of a particle from 
the metastable potential barrier under the influence of a high 
frequency field is shown. The incident radiation can excite 
the tunneling field to higher energy states with higher rate of 
tunneling; however, the probability for being excited to higher 
energy states would decreases. The dotted arrow indicates 
tunneling with dissipation where the field loses some energy 
to the environment. The emerging particle has less energy 
compared to the energy before tunneling and has to travel 
longer distance under the barrier. 


If the condition Sq ^ ai is not satisfied, in addition to 
linear expansion in 5, higher order corrections need to be 
considered. 


In this semi-calssical method the probability of multi¬ 
photon processes in the enhancement of the quantum 
tunneling has been automatically included, [271 [M]. In 
a multi-photon process, the charge trapped in a local 
minima can absorb either one or more than one pho¬ 
ton and tunnels away. Although, absorbing more than 
one photons is less probable; however, upon absorption 
of more photons the particle becomes exited to a higher 
level from which tunneling is facilitated, thus t he pr oba- 
bility of quantum tunneling increases (see Eig. A.9). 
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